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ABSTRACT 

We report on the discovery of four ultra-short period (P^O.18 days) eclipsing M-dwarf 
binaries in the WFCAM Transit Survey. Their orbital periods are significantly shorter 
than of any other known main-sequence binary system, and are all significantly below 
the sharp period cut-off at P ~ 0.22 days as seen in binaries of earlier type stars. The 
shortest-period binary consists of two M4 type stars in a P = 0.112 day orbit. The 
binaries are discovered as part of an extensive search for short-period eclipsing systems 
in over 260,000 stellar lightcurves, including over 10,000 M-dwarfs down to J=18 mag, 
yielding 25 binaries with P ^0.23 days. In a popular paradigm, the evolution of short 
period binaries of cool main-sequence stars is driven by loss of angular momentum 
through magnetised winds. In this scheme, the observed P ~0.22 day period cut-off 
is explained as being due to timescales that are too long for lower-mass binaries to 
decay into tighter orbits. Our discovery of low- mass binaries with significantly shorter 
orbits implies that either these timescales have been overestimated for M-dwarfs, e.g. 
due to a higher effective magnetic activity, or that the mechanism for forming these 
tight M-dwarf binaries is different from that of earlier type main-sequence stars. 
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1 INTRODUCTION 

The period distribution of close binary star systems (e.g. 
Devor 2005; Derekas et al. 2007) contains important infor- 
mation on binary formation and evolutionary processes. Ob- 
servations have revealed that there is a sharp cut-off in the 
period distribution at ^0.22 days (e.g. Rucinski 1992, Nor- 
ton et al. 2011), and very few binaries have thus far been dis- 
covered with significantly shorter periods (GSC 2314-0530 
0.192 days; Dimitrov & Kjurkchieva 2010 and OGLE BW3 
V38 0.198 days; Maceroni k Montalban 2004). Searching 
for binaries beyond this cut-off is interesting because their 

* E-mail: nefs@strw.leidenuniv.nl 



frequency of occurence, and the ratio of detached versus con- 
tact binary systems, provide direct constraints on theories 
that model the formation and migration history of low-mass 
stars, a mass regime that has been fairly poorly charac- 
terised so far. It has even been proposed that very short 
period low-mass binaries could be the progenitors of stel- 
lar mergers which may lead to the observed population of 
very hot Jupiters (Martin et al. 2011), and which could ex- 
plain events such as Nova Sco 2008, where a contact binary 
with an orbital period of 1.4 days merged into a single star 
(Tylenda et al. 2010). 

There are several theories that aim to explain the ob- 
served cut-off and the apparent lack of systems beyond it. In 
the first theory, near-contact binaries are formed from ini- 
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tially well- detached systems that undergo angular momen- 
tum loss (AML) via a magnetised wind on Gyr timescales. 
Evidence for such winds, especially in low-mass, short-period 
binary systems, is provided by observations of extensive cool 
spot coverage, rapid rotation in tidally locked orbits, strong 
Ha emission, flares and high activity rates (e.g. Morales et 
al. 2010, Vida et al. 2009). Stepien (1995; 2006; 2011) es- 
timated the timescale required for the components of a de- 
tached binary to both reach Roche-lobe overflow through 
AML, i.e. to become a contact binary. He proposed that the 
AML timescale is much longer for low-mass systems, such 
that M-dwarf binaries (except in the most extreme mass- 
ratio cases) can not reach Roche-lobe overflow within the 
age of the Universe. The 0.22 day cut-off therefore corre- 
sponds to a lower limit in total binary mass of ~1.O-1.2M0 
for contact systems. However, detached low-mass ultra-short 
period binary systems are also extremely difficult to form 
by this mechanism. In Section 2, we will further explain and 
quantify the AML evolution for M-dwarf binaries. 

Proposing an alternative theory, Jiang et al. (2011) sug- 
gest that an instability in mass transfer, when the primary 
ffUs its Roche lobe, is responsible for the observed short- 
period cut-off. This instability is predicted to occur if the 
primary star has a signiflcant convection zone, implying that 
binaries with a primary mass lower than O.GSM© would 
merge too quickly to form stable contact systems. In this 
scheme, M-dwarf binaries in the contact phase are short- 
lived and are therefore extremely rare. Unlike the AML 
model of Stepien (2011), detached low-mass short-period bi- 
naries are permitted in the mass transfer instability model 
because its timescale for angular momentum loss is signifi- 
cantly shorter. 

However, third-body interactions in the birth environ- 
ment may be responsible for accelerating the orbital evolu- 
tion of low-mass binaries beyond the previous predictions. 
Here, energy is drained from the binary by either ejection 
of the lowest-mass companion in a low-N 'mini-cluster' (e.g. 
Reipurth & Clarke 2001, Goodwin et al. 2004) or through 
the Kozai mechanism (e.g. Eggleton & Kiseleva 2001, Fab- 
ry cky & Tremaine 2007). Interestingly, an adaptive optics 
search for companions around solar-mass contact systems 
by Rucinski et al. (2007) found that the fraction of short- 
period binaries in hierarchical triples is at least 59% (±8%) 
(see also Tokovinin et al. 2006), indicating that a significant 
amount of short-period binaries are found in triple systems. 

With these competing scenarios in mind, it is interest- 
ing to characterise the low-mass M-dwarf population of ul- 
tra short-period binaries. Although M-dwarfs form the most 
common stellar population in our Galaxy (^-^70% by num- 
ber; Henry et al. 1997), their intrinsic faintness is a chal- 
lenge when trying to obtain a sufliciently large M-dwarf 
sample from a magnitude limited optical survey in order to 
study their binary characteristics. In this paper we describe 
a search for the shortest period eclipsing binary systems in 
the WFCAM Transit Survey (WTS) down to J=18. The 
WTS is an infrared photometric monitoring survey running 
on the 3.8m United Kingdom Infrared Telescope (UKIRT) 
with its main goal to flnd planets transiting M-dwarfs. Be- 
cause M-dwarf spectral energy distributions (SEDs) peak in 
the infrared, the WTS is sensitive to redder eclipsing sys- 
tems, with a significant population of binary systems down 
to ~M5. 



In Section 2 we will expand on our motivation of this 
work. In Section 3 we describe our observations and data 
reduction, including the WTS survey itself and low resolu- 
tion spectroscopic follow-up of some of the candidates. In 
Section 4 we use a variability statistic to flnd stars with 
correlated variability and select eclipsing binary candidates 
from the WTS data. We subsequently determine orbital pe- 
riods using box-fltting and Fourier techniques. With simple 
colour cuts we then pre-select binaries with primary mass 
Ml < 0.7 Mq (K5). In Section 5 we obtain estimates of the 
binary effective temperature through broadband SED fltting 
and spectral template matching. By fltting Fourier series to 
the data we identify (semi-) detached and contact systems. 
In Section 6, we compare our results with expectations from 
the evolution scenarios of low-mass short period binaries. 



2 MOTIVATION 

2.1 The AML timescale argument 

Single stars with masses ^ 1.5^4© have an outer convec- 
tive layer and can be chromospherically very active. This 
activity is expected to cause the loss of angular momentum 
through a magnetised wind of material that is forced to co- 
rotate with the star. The rotation rate of a star is found to 
correlate with the activity. Young active stars rotate rapidly 
but their rotation slows down with time leading to a de- 
crease in their activity (e.g. Irwin et al. 2011). In a close 
binary, with synchronised orbits and spins, the angular mo- 
mentum loss due to the spin-down of the individual stars 
causes a decrease of the binary orbital angular momentum, 
which tightens the orbit, but this then spins up the stars. A 
popular scenario is that when given enough time, the loss of 
angular momentum will bring the two stars into contact. 

Stepien (1995, 2006, 2011) estimated timescales re- 
quired for a detached binary to reach Roche-lobe overflow 
through AML and to evolve into a contact system, conclud- 
ing that for a solar-type, equal-mass binary with a starting 
separation that corresponds to a period Po=2.0 days this 
timescale is '-^6.5Gyr. He proposed that the AML timescale 
is signiflcantly longer for decreasing mass systems and that 
therefore Roche-lobe overflow will not occur within the Hub- 
ble time for binaries with an initial primary mass of less than 
~O.7M0, again assuming an initial orbital period of 2.0 days 
(an assumption which is further explained below). The 0.22 
day period cut would therefore correspond to a lower limit on 
the total mass of contact systems of 1.O-1.2M0 (i.e. a binary 
consisting of two MO dwarfs). This suggests that evolution 
timescales are too long for M-dwarf binaries to decay into 
tight orbits within the age of the Universe. 

One solution to this is to choose a shorter starting pe- 
riod {Po). To calculate the orbital evolution for M-dwarfs, 
and to investigate this evolution for different starting pe- 
riods and mass ratios, we equate the change in the total 
binary orbital angular momentum (Horb), which scales with 
the orbital frequency cj as dHorb/dt ~ cj^^^^cj, to the change 
in the sum of spin angular momenta (Hspin) for the binary 
components, dHspin/dt ^ cj, following the derivation in Sec- 
tion 3 of Stepien (1995). We assume that mass- loss is small 
over the Hubble time for M-dwarfs and that the orbit re- 
mains circular and synchronised throughout binary evolu- 
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Figure 1. Panel a): The orbital evolution of M-dwarf binaries, assuming Po=2.0 days and q=1.0, plotted for decreasing primary stellar 
mass from 0.7 to O.2A40 in steps of O.IAI0 (black solid curve). The dashed grey curves are models with mass ratio q=0.3. The models 
are cut off at the point of Roche lobe overflow, which depends on total mass and the mass-ratio. The horizontal line indicates the 0.22 
day period cut-off, whereas the vertical line shows the Hubble time. Binaries with mass less then ~O.65A^0 could not have reached 
contact. Panel b): same models, but for Po=1.0 day. Panels c) through e): the timescales to reach contact as a function of primary 
mass for Po=2.0, 1.5 and 1.0 days and mass ratios q=1.0, 0.7, and 0.3. 
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tion. The change in the orbital period as function of time 
can then be written as: 
r7 P 

^ = -A{r^,Mi,r2,M2)P-'^\ (1) 

where A is constant with time and ri,2 and Mi, 2 are the 
radii and masses of the primary and secondary binary com- 
ponents. This equation has the solution: 

Pit) = ip^^' - ^Atr^\ (2) 

where Po is the starting period at formation (in days) and 
P(t) is the observed orbital period at time t. 

In panel a) of Figure [1] we show the orbital evolution 
of M-dwarfs, assuming Po=2.0 days and q(=A^2/A^i) = l-0, 
plotted for decreasing primary stellar mass from 0.7 to 
0.2A4q in steps of O.IMq. In the same panel, we show the 
evolution for M-dwarfs with q=0.3 (grey dashed curves). Bi- 
nary evolution is significantly faster for low-mass ratio sys- 
tems. Clearly, unless a low mass-ratio is assumed, there is 
insufficient AML over the Hubble time for any M-dwarf bi- 
nary to reach contact. In panel b) of Figure [T] we plot the 
same binary evolution for M-dwarfs but with Po=l-0 day, 
showing that a shorter starting period results in faster evo- 
lution. 

In panels c) through e) of Figure [T] we show the time 
required to reach contact as a function of primary mass as- 
suming Po=2.0, 1.5, and 1.0 days for mass ratios q=1.0, 
0.7, and 0.3. Clearly, in the Po=l-0 day model, evolution 
is much faster and q=l systems with primaries more mas- 
sive than ^0.35Mq could be brought into contact within 
the Hubble time. This means that the upper limit on the 
total mass of contact systems would shift to masses cor- 
responding to later-type M-dwarfs, assuming that the for- 
mation mechanism for such binaries is the same as that of 
earlier-type stars. Also, even for a O.S+O.SA^© binary the 
evolution timescale is ^8 Gyrs in the fastest model, which 
would mean that such binaries would be some of the oldest 
objects in our current universe. 

We have calculated what Po is minimally required to 
bring a given M-dwarf binary into a given orbital period 
today, within a given time- frame. We show this as a func- 
tion of primary mass, derived using Equation 2, in Figure 
[2] for equal-mass binaries, assuming observed orbits of 0.10, 
0.15, 0.20 and 0.25 days. For example, a minimum starting 
period of Po = 1.0 day is required for a O.3+O.3A^0 binary 
to reach P=0.2 days within the Hubble time. For a typical 
5.0 Gyr thin disk binary, the required Po is 0.4 days for the 
same system. Clearly, to constrain AML theory, it is vital 
to assess from observations what is the actual frequency of 
(ultra) short-period M-dwarf binaries. 

Although shorter starting periods could at least partly 
explain the existance of ultra-short period M-dwarf systems 
within the AML framework, binaries with these short start- 
ing periods should then be found abundantly in young clus- 
ter environments. However, pre- main sequence binary pro- 
genitors with Po ^2.5 days appear to be very rare, with 
4 currently observed. HD15555 (a G5+K1 in a 1.7 day 
orbit), MML53 (a G2+K2 in a 2.1 day orbit), V4046Sgr 
(K5+K7;2.42 days) and 155913-2233 (K5+K5;2.42 days) are 
the only known (Melo et al. 2000, Hebb et al. 2010). Also, 
contact binaries of any spectral type are extremely rare in 
open clusters younger than ^4 Gyr, only TXCnC (0.38 day 
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Figure 2. The starting period Pq required to bring an equal-mass 
binary of given primary mass into a currently observed orbit of 
0.1,0.15,0.2, and 0.25 days (the four lines, observed period in- 
creasing from bottom to top) for t=13.8 Gyr (upper set of four 
lines) and t=5.0 Gyr (lower set of four lines). The thick black 
solid curve indicates the estimated minimum possible separation 
period at birth for a binary using 1.0 Myr Baraffe models. 

orbit; Zhang et al. 2009) is linked to the intermediate-age 
open cluster Preasepe (^600 Myr). Cargile et al. (2008) re- 
port the discovery of a 4.7 day eclipsing pre-main sequence 
M-dwarf with masses 0.39 and 0.40^4©, which suggests that 
shorter-period pre-main sequence M-dwarfs are not subject 
to observational bias. 

The minimum possible orbital separation of a close bi- 
nary during the fragmentation of a protostellar cloud is an 
important constraint on the assumed Po as it prevents ex- 
plaining ultra-short period binaries as simply forming very 
close together. This is because T-Tauri stars of subsolar mass 
have radii of ^2-3R0 (Baraffe et al. 1998). A pile- up of de- 
tached binaries with periods of 2-3 days is also expected 
because the Kozai mechanism in a triple system can effec- 
tively produce such binaries within 50 Myrs (e.g. Fabrycky 
and Tremaine 2007), but is less efficient for shorter periods. 
We approximately estimate the shortest possible Po assum- 
ing pre-main sequence stellar radii using Baraffe (1998) 1.0 
Myr models (the thick black solid curve in Figure [2|). Note 
that for the lowest mass M-dwarfs no Po can yield very short- 
period binaries. This makes it difficult to justify any starting 
period of Po ^ 1.5 days. 

2.2 Periods at contact 

We estimated the typical orbital periods we would expect 
for M-dwarf contact binary systems in order to focus our 
search for these systems. To do this, we determined their 
critical Roche lobe size ri,2 (e.g. Eggleton 1983). For exam- 
ple, for an equal-mass binary this is approximately given by 
r\ — T2 ^ 0.379a, where a is the binary separation. Using 
the main-sequence mass-radius relation for low mass stars 
at IGyr of age (Baraffe et al. 1998) and Keplers law we 
calculate corresponding critical Roche lobe periods Pnoche- 
For illustration this Pnoche is plotted as a function of pri- 
mary mass in Figure [3] (solid curve) . For a binary of two 
MO stars (assuming a primary mass of O.6M0; Baraffe & 
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Figure 3. The estimated binary orbital period (Roche lobe pe- 
riod) for low-mass q=l contact binary systems, plotted as func- 
tion of primary mass for three different models, all assuming the 
Baraffe et al. (1998) IGyr mass-radius relation: i) a spherical 
model without distorted radius (solid line); ii) a fully radiative 
n=3 polytrope model (dashed line;no distortion) and iii) a fully 
convective stellar atmosphere with distortion due to tidal effects 
and rotation (Knigge et al. 2011; dotted line). The thick vertical 
line indicates the full convection limit (^O.SSA^©), whereas we 
also indicate approximate KO, M0,M3 and M5 type stars. Also 
shown are the 0.22 day period cut (horizontal solid line) and the 
five shortest period (P ^0.21 days) eclipsing binaries from the 
WTS red binary sample (presented in Section 6). 



Chabrier 1996), PRoche is -0.19-0.21 days, for a twin M5 
system —0.05-0.07 days and —0.35 days for a solar type bi- 
nary. The Roche period is also a function of q; for a 0.6Mq 
system with q=0.8, PRoche is —0.13 days and —0.07 days 
for a q=0.5 system. This indicates that binaries with lower 
q have even shorter PRoche • The dotted curve in Figure [3] 
is computed from the expression for the radius of a tidally 
and rotationally deformed Roche-lobe filling star (Knigge 
et al. 2011, Section 5.2.2). For such a star PRoche is higher 
because the equilibrium radius of each star is now affected 
by the deformation. Note that there is still considerable un- 
certainty in the mass, radius, Tgff, spectral type calibration 
(e.g. Baraffe & Chabrier 1996), and magnetic activity may 
be a contributing factor to radius inflation in a tight low 
mass binary system (e.g. Rozyczka et al. 2009), which in 
turn increases the length of PRoche- 

Therefore, to critically test the predictions of AML and 
other paradigms, it is essential to determine from observa- 
tions the abundance of M-dwarf binaries that are near con- 
tact, i.e. binary stars which are expected to have ultra-short 
orbital periods on the order of only a few hours. Because 
of the intrinsic faintness of M-dwarfs and the strong line 
broadening induced by rapid stellar rotation it is difficult 
to find such binaries using data from current radial velocity 
surveys. We therefore present the results of a campaign util- 
ising photometric measurements obtained by the WFCAM 
Transit Survey, which has both the cadence and sensitivity 
to detect and characterise ultra-short period binaries in the 
uncharted M-dwarf regime. 



Name 


(RA) 


(DEC) 


#stars 


#epochs 


#cand. 


03hr field 




+ 39^^14"^ 


10827(36306) 


392 


74 


07hr field 






16623(56070) 


626 


140 


17hr field 


17^14"^ 


+ 03^^44"^ 


9621(39879) 


709 


68 


19hr field 


19^35"^ 


+ 36^^29"^ 


34452(130320) 


1154 


375 



Table 1. The main properties of the four WTS survey fields. In- 
dicated are the approximate centres of the fields (right ascension 
and declination), the total number of stars with J ^16 (J ^18 in 
brackets), the number of epochs as of September 7th 2011, and 
the number of binary candidates per field. 



3 OBSERVATIONS AND DATA REDUCTION 

3.1 WTS J band time-series photometry 

The WFCAM Transit Survey (WTS) is an ongoing photo- 
metric monitoring campaign operating on the 3.8m United 
Kingdom Infrared Telescope (UKIRT) at Mauna Kea, 
Hawaii, since August 2007. The Survey uses the Wide- 
Field Camera (WFCAM), which has four 2048x2048 18/im 
HgCdTe Rockwell Hawaii-II, infrared imaging arrays which 
each cover 13.65' x 13.65' (0.4" /pixel), and are separated by 
94% of a chip widttQ (Casali et al. 2007). Observations for 
the WTS are obtained in the J-band (1.25/im), near the 
peak of the spectral energy distribution (SED) of a typical 
M-dwarf. 

The WTS is primarily designed to find planets transit- 
ing red dwarf stars and to detect low mass eclipsing binaries, 
by observing ^-^6000 early to mid M-dwarfs with J ^16 (in 
this paper we consider stars down to J=18). The WTS ob- 
serving strategy uses the queue-schedule mode of UKIRT 
and can operate in mediocre seeing and thin cloud cover. 
Four target fields, 1.5 square degree each and passing within 
15 degrees of zenith, were selected to give year-round visibil- 
ity. The fields were selected as regions of the sky where the 
ratio of dwarfs to giants was maximized, reddening relatively 
low (E(B-V) between 0.057 and 0.234) and overcrowding re- 
duced by observing close to but outside the galactic plane 
{b > 5°). In Table [1] we summarise the main properties of 
the four survey fields (the 03hr field, 07hr field, 17hr field 
and 19hr field). A field is observed by dithering the four 
detectors of WFCAM (a paw print) through 8 pointings, 
where each pointing is labeled with a single letter from a to 
h (e.g. '19a'), which are then tiled to give uniform coverage 
across the field. In itself, each single pointing consists of a 
9-point jitter pattern with 10 second exposures at each jit- 
ter position. In this way, the near-infrared lightcurves have 
an average cadence (including overheads) of 15 min. Typi- 
cally, WTS observations are taken only at the beginning of 
a night, just after twilight in >1" seeing. 

Data reduction of the raw 2D J-band images is per- 
formed using the Cambridge Astronomical Survey Unit 
(CASU) pipeline, which is based on the INT wide-field sur- 
vey pipeline (Irwin & Lewis 2001). Astrometric and photo- 
metric calibration is obtained using 2MASS. Source detec- 
tion is done on stacks of the 20 best exposures (the 'master 
frame') and aperture photometry is performed as described 
in Irwin et al. (2007). The brightest stars (saturation oc- 
curs around J= 12- 12.5) have a precision of ^-^3 mmag per 
data point, whereas 1% photometry is still achieved for J=16 
(-7% for J=18). 



http : / / casu . ast . cam .ac.uk/ surveys-pro j ects / wfcam 
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The short-period ecUpsing binary candidates described 
in this paper were obtained from WTS hghtcurves reduced 
on 7 September 2011 of all four survey fields. The 19hr field 
had the most extensive coverage (^1100 data points). For 
each field, we also obtained single deep exposures in WF- 
CAM ZYJHK. These wih be used together with griz pho- 
tometry from the Sloan Digitized Sky Survey (SDSS) where 
available to create SEDs and estimate the effective temper- 
ature of the candidates. 



3.2 INT broad-band photometry 

The 17hr field in the WTS is not covered by SDSS, therefore 
we obtained optical photometry of this field to use in our 
SED fitting and effective temperature estimates. We used 
the Wide Field Camera (WFC) on the 2.5m Isaac Newton 
Telescope (INT) in La Palma on the nights of 8,10 July 
2009 and 26 July 2010 to obtain g, r and i photometry of 
this field. The smaller field-of-view of WFC compared to 
WFC AM required 14 pointings to cover the entire field. We 
used the following exposure times for each pointing: g 2x 250 
seconds, r 500 seconds, and i 100 seconds, in order to match 
the typical depth achieved by SDSS. The data were reduced 
using custom-made IDL procedures to bias-subtract, flat- 
field (and de- fringe in the case of the i-band, where the 
fringe-frame was created from a number of deep exposures 
of the moonless sky). We used the IDL procedures find.pro 
to identify the sources and aper.pro to extract the photome- 
try, calibrating to SDSS magnitudes using INT observations 
of the 19hr field in the same filters. 



3.3 Low resolution spectroscopy 

As an independent check on the effective temperatures de- 
rived for our binary candidates in Section 5.2, from broad- 
band colour fitting of their SEDs, we obtained low resolution 
spectral observations of three of our targets. Optical spectra 
of 19h-3-14922 (P=0.1798 days, (r-i) = 1.65) and 19g-3-06701 
(P=0.2090 days, (r-i) = 1.01), were taken with the 4.2m WHT 
on the night of August 17 2010. We used the red arm and 
the R158 grism of the ISIS dual beam spectrograph, giving a 
dispersion of ^1.8A/pixel using a slit width of 1.0" matching 
the seeing. In addition, a spectrum of 07e-2-03887 (P=0.2734 
days, (r-i) = 1.0) was obtained with ACAM (^3.3A/pixel) 
on the night of 08 September 2011. The spectra were re- 
duced using a combination of IRAF (the apall package) and 
custom IDL routines. In IDL, the spectra were trimmed to 
encompass the length of the slit, bias-subtracted, flat field 
corrected and median- filtered to remove cosmic rays. Back- 
ground subtraction and optimal extraction of the ID spec- 
tra was performed in IRAF. We used frames taken with the 
CuNe+CuAr standard lamp to determine the wavelength 
solution. For flux-calibration, we obtained reference spectra 
of standard stars. 



4 SAMPLE SELECTION 
4.1 Variability statistic 

The WTS 3.0 release consists of 262, 000 stellar light curves 
with J ^ 18. First, we select only those light curves that have 



Variability index 



14 16 
Median J band magnitude 



Figure 4. The Stetson Jg variability statistic for one detector 
of the 19h WTS paw, plotted as a function of WFCAM J band 
magnitude. We use a cut Js > 0.2 to select candidate eclipsing 
binaries. At the brightest magnitudes, outliers are dominated by 
systematic effects. These are removed through visual inspection. 



less than 20% bad datapoints and clip the bad data before 
further processing. We define as bad data a J-band measure- 
ment that either has infinite value (saturation, near-edge 
detections) or a photometric error larger than 1.0 magni- 
tudes. We subsequently use the Stetson Js'-statistic (Stetson 
1996) as a quant at ive estimator for photometric variability, 
analogous to Pepper & Burke (2006). The Stetson statistic 
weights photometric observations that are correlated in time 
using subsets of data that are separated by less than 0.03 
days. For a nonvariable star showing only random (Gaus- 
sian) noise, this statistic will be around 0, whereas it will be 
positive for stars with correlated (physical) variability. As 
an example in Figure |4] we show the Js variability index for 
one WFCAM detector in WTS paw 19h. To select candidate 
variable stars, we apply a conservative cut of Js,cut >0.2. We 
determine this cut by examining the length of our candidate 
lists for random WTS paws for values of Js,cut from to 0.5. 
We subsequently checked each candidate lightcurve and ob- 
tain a cutout thumbnail image from the stacked WFCAM 
J band master frame to eliminate the most obvious false- 
positives through visual inspection. These false-positives are 
caused by blending effects from nearby stars, extended wings 
and bleed trails from bright stars or detector defects. This 
step eliminates 80-90% of the candidates and typically 20 to 
40 objects per individual detector pass the inspection stage. 
In total, we find 656 variables with J ^18 in this way. 375 
candidates are from the 19hr field (which contributes 49.6% 
of the total number of hghtcurves with J ^18), 68 from the 
17hr field (15.2%), 140 from the 07hr field (21.4%) and fi- 
nally 74 candidates are from the 03hr field (13.8%). In Figure 
Owe show a plot of the candidate fraction as a function of J 
band magnitude. At J=16, roughly 0.5% of the high quality 
hghtcurves is classified as a candidate, but for fainter mag- 
nitudes this fraction quickly decreases, reflecting the poorer 
photometry. For this reason, we have not investigated any 
hghtcurves with J ^18. 
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Figure 5. The ratio of binary star candidates over the number 
of available good quality WFCAM lightcurves as a function of 
WFCAM J magnitude. Each black dot indicates this fraction de- 
termined in 0.1 magnitude bins. The plot is truncated at J=12 
because of saturation effects for brighter stars. 

4.2 Orbital period determination 

To determine periods for all binary candidates we adopt a 
simple box- fitting algorithm (Kovacs et al. 2002), such as 
used by Collier- Cameron et al. (2006). We limit our search 
to P ^0.1 days because our frequency spectra are very noisy 
at the shortest periods, which we suspect is due to strong 
aliasing caused by the peculiar form of the WTS window 
function. Note that shorter period systems could still be 
detected as an alias of the real period. We phase- fold the 
lightcurves over 20,000 periods between 0.1 and 1.0 days and 
for every period we vary the primary eclipse epoch To and 
box width, and calculated the error- weighted signal com- 
pared to white noise (S/N). We then refined our analysis by 
selecting the period and epoch with the highest significance 
and re-iterate the box-fitting procedure with a smaller range 
of parameters, centered around the frequency peak. In this 
way, we generally achieve uncertainties in period of better 
than ^ 10~^ days and better than ^ 10~^ days on Tq. In 
this procedure we exclude periods within 10 minutes of 0.33, 
0.5 and 1.0 days. In the case the algorithm confuses the sec- 
ondary with the primary eclipse, this is picked up through 
visual inspection. In these cases the lightcurves are folded on 
several period aliases. As a crosscheck on the performance 
of the box- fitting algorithm, we also compute the Lomb- 
Scargle power spectrum (Scargle 1982) between 1 hr and 2 
days (with a period sampling of 0.0001 days) for each candi- 
date. Since the Lomb- Scargle method uses Fourier analysis, 
it is mostly sensitive to W Uma variables, pulsators and 
rotators, and less sensitive to the narrow eclipses that are 
expected for (semi-) detached systems. In the end we adopt 
the method that provides the lowest rms in the phase-folded 
light curve. 

4.3 Selection of the final sample 

We applied a first-pass colour cut using SDSS photometry 
(Kowalski et al. 2009) to select late K and M-dwarf short- 
period eclipsing binaries, which is verified using the fitting 




Figure 6. SDSS colour-colour diagram for objects in the 19hr 
field, showing (r-i) versus (i-z). The light grey solid lines indicate 
the median average M-dwarf spectral types from Kowalski et al. 
(2009). Note that these averages generally have la uncertainty 
of 0.05 to 0.15 in (r-i) and (i-z). Black dots are objects in the 
WTS 19hr field (for plotting purposes we only show one field), 
which have SDSS DR7 photometry. Dark grey squares indicate 
our 'blue' short period EB sample (P ^0.23 days), dark grey 
filled dots indicate our red sample, where the light grey stars in 
the figure are red sample binaries (with 0.23 days). 



of the broadband SEDs in Section 5.2. SDSS photometry is 
available for three WTS survey fields, the 03hr field, 07hr 
field and 19hr field. For the 17hr field, which does not have 
SDSS coverage, we use gri photometry obtained on the INT. 
We made cuts at (r-i) ^0.50 and (i-z) ^0.25 and selected or- 
bital periods with P ^0.3 days. This defines our 'red' sam- 
ple, the principal focus of this paper. In addition we use the 
well-calibrated SDSS colour selection criteria from Sesar et 
al. (2007) to exclude RR Lyrae. As there may be overlap in 
estimated effective temperature between late K-type stars 
and early M stars we also kept binaries which have P ^0.23 
days (~20,000 seconds, analogous to the sample of Norton 
et al. 2011) that are bluer than our cut, but are excluded as 
pulsators using the Sesar et al. (2007) colour criteria. This is 
our 'blue' sample. We have removed objects from our red and 
blue samples that have sinusoidal phase-folded light curves 
with a large scatter in their amplitude, which would indicate 
variability caused by star spots (the amplitude depends on 
the spot-size which is expected to vary on timescales shorter 
than the 4 year observing span of the WTS). Note that this 
could remove a small fraction of genuine grazing equal-mass 
contact systems at our faintest magnitudes where our pho- 
tometry is less precise. 

In Figure [6] we show our red and blue samples in a (r-i) 
versus (i-z) SDSS colour-colour diagram. The light grey lines 
indicate the median M dwarf stellar types, as presented in 
Kowalski et al. (2009). The colours of these spectral types 
generally have a la uncertainty of 0.05-0.15. Dark grey filled 
dots and light grey stars indicate objects from our red sam- 
ple, whereas dark grey filled squares represent the blue sam- 
ple. The lightcurves and the optical colours of all of our 
candidates do not show features which would suggest they 
are M-dwarf white dwarf systems (Covey et al. 2005). 
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Figure 7. Top panel: Classification of our binary sample in the 0,2, 0-4 Fourier coefficient plane, where the numbers identify the 
candidates of Table 2, and filled squares refer to the red sample whereas asterisks indicate the blue sample. The lower solid line denotes 
the envelope for contact, which is relating the two coefficients through a4 = a2 (0.125 + a2). Sources below the envelope are classified 
as contact binaries. 'SD' indicates semi-detached binaries. Lower panel: two representative binaries classified as detached and contact 
systems. The grey solid line is the best- fit Fourier decomposition (Section 5). 
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4-3.1 Reddening estimates 

We have investigated the reddening effect for the 19hr 
field using a model for interstellar extinction presented by 
Amores & Lepine (2005). At J ^ 17, assuming no redden- 
ing, the WTS is distance- limited to ^^1800 pc for the earliest 
M-dwarfs. At 1.8 kpc, the reddening model gives an extinc- 
tion of Av=0.37 mag (E(B-V))=0.12 mag) in the direction 
of the 19hr field, and using relations presented in Table 6 of 
Schlegel et al. (1998), we converted between filters to find 
Ag=0A3 mag, Ak=0M2 mag, E(r-i)=0.076, E(i-z)=0.069 
and E(J-H) =0.037. The reddening effect along the line-of- 
sight to the field thus appears to be small, and similar results 
are obtained for the other three WTS fields. This makes it 
unlikely that any of our short-period red binary candidates 
are highly reddened pulsators. Also note that the pulsation 
periods of giant stars are inversely proportional to the square 
root of the mean stellar density, meaning that the low densi- 
ties of giant stars result in significantly longer pulsation pe- 
riods than we detect in our candidates. Rodriguez-Lopez et 
al. (2012) suggested that M-dwarfs may be pulsating, how- 
ever it is not clear if such M-dwarf pulsations are stable and 
whether they could be detected in M-dwarfs with Gyr ages. 
They have not been confirmed observationally to date. Note 
that any pulsation would also have to be stable on the order 
of four years (the span of the WTS observing campaign) . 
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Figure 8. Two low resolution WHT ISIS red arm spectra and one 
ACAM spectrum of our short period binaries. The spectra are flux 
calibrated and normalised. Top panel: 19h-3-14922, a detached 
system of ~M3 and M5 stars. Middle panel: 07e-2-03887, a 
(semi-)detached system of M2 stars. Lower panel: 19g-3-06701 
a contact system of ~M0 and Ml stars. The grey curves are best- 
fit NextGen model spectra for solar metallicity, log(g)=5.0 and 
age 1.0 Gyrs. The dotted vertical line indicates the position of the 
Hq; line at 6563 A, which is related to emission from the stellar 
chromosperes induced by spun up rotation in the binary. 



5 CHARACTERISATION OF THE ECLIPSING 
BINARY SYSTEMS 

5.1 Binary classification 

For a crude initial classification of a binary in the three cat- 
egories detached, semi-detached or (over) contact, we follow 
Rucinski (1993) and decompose the phased lightcurve m($) 
into Fourier components with coefficients ai and bi: 

10 

m($) = mo + ^[tti cos (27ri$) + bisin{27Ti^)], (3) 

where mo is the mean lightcurve magnitude. The number 
of Fourier components is limited to 10, because for higher 
i the quality of the fit does not improve for the typical un- 
certainties reached. For the fitting of the Fourier coefficients 
we use the IDL routines fourfit and fourfunc from Marc 
Buie's IDL archival. Based on a grid of coefficients derived 
from fits to model eclipsing binary lightcurves the a2 and 
a4 coefficients allow a distinction to be made between the 
three categories by plotting the sources in the (a2, a a) plane 
(Figure [T)). The coefficient a2 is related to the depth of the 
eclipses, which depends on the inclination of the orbit and 
the mass ratio. For example, we find a2 ^^0.036 for 19b-3- 
06008, a relatively shallow semi-detached binary, whereas 
deep eclipse systems such as 07c-4-05645,19e-4-00861 and 
19c- 2-08140 have a2 - 0.13-0.16. The coefficients ai and as 
are related to the ratio of primary versus secondary eclipse 
depths and bi to the difference in the lightcurve maxima at 
phase 0.25 and 0.75 (star spots or mass transfer can bias 
this value). Pulsating variables can be distinguished from 

Available at |http: / / www.bQulder.swri.edu/ ^buie/idlTl 



genuine eclipsing systems, in addition to the colour cuts pre- 
viously described, by using the 62 and 64 coefficients (e.g. 
Pojmanski 2002). Our candidates have 62 and 64 consistent 
with non-pulsating stars. 

5.2 The effective temperature 

We estimate the effective temperature Teff of the sources in 
our sample using an 8-band weighted SED fit to combined 
SDSS griz and WFCAM YJHK photometry. We fit the 
observed SED to a grid regularly spaced by 50 K in Teff with 
theoretical broad-band colours based on Baraffe et al. (1998) 
IGyr, log(g)=5.0, solar metallicity models. In this process, 
we re-iterate the fitting procedure after removing the band 
that gives the largest outlier from the initial fit. For three 
sources we obtain an independent check on our results using 
low resolution spectra from the WHT (Section 2.2). For the 
spectroscopy data we fit NexTGen model spectra (Allard 
et al. 1997) with steps in effective temperature ATeff=100K, 
solar metallicity and surface gravity log (g) =5.0 to our ob- 
served spectra resulting in uncertainties in spectral type of 
about one subclass (corresponding to ^100-150K). We mask 
out the strong telluric ozone bandhead around 7600 A and 
the He emission line at 6563 A. 

In Figure [8] we show our three low resolution spectra 
(solid black lines) with overplotted best-fit model spectra 
(grey lines). In this figure we indicate approximate spectral 
types for the binary components by crude estimates based 
on the simplified assumption that the stars radiate as black- 
bodies. We use the relative depths of the primary and sec- 
ondary eclipses to obtain the ratio of surface brightness and 
together with the ratio of radii, which we estimate from the 
square of the primary and secondary eclipse depths, we solve 
for the individual temperatures Ti,2, analogous to Goughlin 
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et al. (2011). For the conversion to spectral type we use the 
tabulated model spectral types as a function of Tgff from 
Baraffe & Chabrier (1996). Note that the models consis- 
tently underpredict the flux around 6900 A, which may be 
related to missing opacities of TiO (Titanium-Oxide) in the 
code we use. We find that the Teff derived from the SED fit- 
ting is generally lower than the corresponding temperature 
from the spectral template matching, typically by one to two 
subclasses, corresponding to 100-300 K, which we attribute 
to missing opacities in the optical. However, given that we 
estimate typical errors of 200-300K on the SED result, the 
two methods are still in agreement. In Table 2 we show Teff 
for each binary candidate. 

All three binaries show significant emission, which 
indicates chromospheric activity (the absense of the Lithium 
absorption line appears to rule out accretion from young 
stars). Using the IDL routine feature we estimate the equiv- 
alent width (EW) of the Ha- The calculated EW values ( 
-7.9 A for 19h-3-14922, -3.8 A for 07e-2-03887 and -5.7 A 
for 19g-3-06701) are consistent with or slightly higher than 
those corresponding to active single stars with similar spec- 
tral types (Riaz et al. 2006). The three systems have equiv- 
alent widths between 30-70% of the empirical accretor/non- 
accretor division for their estimated spectral types (Bar- 
rado y Navascues & Martin 2003). The presence of the Na 
(Sodium) absorption doublet around 8200 A and the absence 
of a strong Ca (Calcium) triplet in the near-IR indicate that 
the three candidates are high gravity objects, i.e. M-dwarfs 
and not giants. 



6 RESULTS & DISCUSSION 

The final sample consists of 25 eclipsing binary systems with 
orbital periods shorter than 0.23 days (20,000 seconds), of 
which 9 have optical to near-infrared colours or optical spec- 
tra consistent with M-dwarfs. Four of these 9 have periods 
shorter than any known main-sequence binary. In Figure 
[9] we show the WFCAM J band lightcurves of these four 
binaries phase-folded to the best-fit epoch and period de- 
termined from the box- fitting algorithm. Five binaries from 
the red sample with periods shorter than 0.23 days have 
Fourier coefficients (a2,a4) typical of a contact system ('C'), 
two binaries are identified as semi-detached systems ('SD'), 
whereas two other systems are detached ('D). From the blue 
sample, 9 out of 14 binaries 64%) are identified as contact 
systems. We have summarised the binary characteristics for 
both the red sample and the blue sample in Table [2l in- 
cluding the Fourier classification. For all of the 25 binaries, 
lightcurves are shown and coordinates are given (in Table 3) 
in the Appendix. 

The shortest period system, 19b-3-06008, is a binary 
consisting of two --M4 stars (Teff --3150 K) with P--0.1121 
days. This binary appears to be close to contact. Assum- 
ing orbital synchronisation, the estimated masses and radii 
for this binary correspond to a fast rotational velocity be- 
tween 15-25% of the breakup speed (Herbst et al. 2001), but 
slow enough to keep the system stable. A second system, 
07g-3-05744, has a period of ^0.151 days, shows slightly 
unequal eclipses and a Teff of 3400 K, suggesting this is 
a binary with spectral types of ^M2-3. This system is 
identified as a contact system by the Fourier decomposi- 



tion. Two other systems, 17d-3-02440 (-Ml. 5) and 19h-3- 
14922 (WHT ISIS spectrum indicating -M3 and M5 stars) 
with periods —0.1572 and —0.1798 days, are detached bina- 
ries with unequal ratio of primary versus secondary eclipse 
depth. Figure [3] illustrates why these binaries are classified as 
detached systems despite their ultra-short periods, whereas 
07d-2-02291 (a blue sample K type binary with P-0.2013 
days) is a contact system; the critical Roche lobe period 
for contact depends on primary mass. Nevertheless, signifi- 
cant out-of-eclipse variability in the lightcurves of the four 
ultra-short period binaries indicates there is significant tidal 
interaction between the primary and secondary in every sys- 
tem. 



6.1 Comparison with previous studies 

This study probes for the first time the low- mass regime 
of ultra-short period eclipsing binaries in an extensive way. 
Recently, Norton et al. (2011) presented 43 eclipsing bi- 
nary candidates identified in the SuperWASP Survey with 
periods shorter than 0.23 days, but only three objects in 
their sample have (J-K)^0.8 (—MO spectral type). Their 
shortest period system is GSC2314-0530 (P-0.1926 days), 
which was modelled by Dimitrov & Kjurkchieva (2010) as a 
semi-detached M-dwarf binary with components of O.51M0 
and O.26M0 (corresponding to spectral types of —MO and 
M4). Two other binaries in their sample have P ^0.21 
days (0.20908 and 0.20964 days). Maceroni & Montalban 
(2004) found an almost twin M3 system (OGLE BW3 V38; 
semi-detached, P=0.1984 days) with strong chromospheric 
emission and masses of 0.44 and O.41M0. Other surveys 
found only binaries with P>0.21 days. Miller et al. (2010) 
found seven candidate contact binaries with periods between 
0.2109 and 0.23 days in a 0.25 square degree Galactic plane 
survey with the ESO 2.2m telescope, and Weldrake et al. 
(2004) found a P=0.2155 day system in the globular cluster 
47 Tuc, whereas Pribulla et al. (2009) report a P=0.211249 
day system. The left panel of Figure [TOl shows a histogram 
of all literature binaries plus our new systems with P^O.23 
days (where grey and black filled bars indicate those from 
the WTS). The WTS sample follows the literature distri- 
bution, but suggests that rather than a sharpe cut-off at 
0.22 days, the population extends into a tail of significantly 
shorter periods for lower-mass systems. 



6.2 The period-colour relation 

Our observations suggest a correlation between spectral type 
and the shortest possible binary period. This is in line with 
the period-colour relation as found by Eggen (1967) for con- 
tact binaries, implying that bluer systems have longer or- 
bital periods. This is because for main-sequence stars the 
effective temperature (colour) is linked with stellar radius, 
and as a consequence orbital period for a contact binary. 
The period-colour relation can be used to identify genuine 
contact systems and weed out pulsating stars with compa- 
rable periods (i.e. 6 Set and 7Dor) and has been used for 
distance estimates (e.g. Rucinski 2007). 

Recently the same relation was presented for infrared 
colours by Deb & Singh (2010) using 2MASS (J-K) pho- 
tometry, which is relatively free from extinction effects. We 
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Figure 9. The (median subtracted) WFCAM J band lightcurves of the four candidate ultra-short period echpsing M-dwarfs binaries, 
phase-folded to the best-fit primary eclipse epoch and orbital period. The boxes in the lower right corner of each plot show a three colour 
SDSS composite image of each source. For 17d-3-02440, which does not have SDSS coverage, we show the WFCAM J band image. A 
colour version of the thumbnails is available in the online paper. 



have collected (J-K) colours from the literature of P ^ 0.23 
days binary systems and combined them with our WTS re- 
sults in a period-near-infrared colour diagram in Figure [TOl 
(right panel). In this figure the stars indicate WTS objects, 
whereas the open boxes are literature sources. We have con- 
verted the WFCAM infrared photometry to 2MASS colours 
using the equations in Hodgkin et al. (2009). We extrapo- 
lated the best-fit powerlaw period-colour relation of Deb & 
Singh (black solid line) to the ultrashort periods of our WTS 
binaries. The dotted curves show the 1 and 3a boundaries 
of this power-law index. There is significant scatter in the 
period versus colour sample, yet, ^65% of all systems with 
P ^0.23 days fall within the la limits of the correlation, in- 
dicating that the majority of our systems are consistent with 
(near-) contact binaries rather than pulsators. However, we 
believe this scatter is intrinsic (the 2MASS J and K bands 
are obtained at the same epoch) and could be due to vary- 
ing degrees of contact or age differences within the sample. 
Our shortest period binary system (19b-3-06008, P^O.1121 
days) is an outlier, indicating that a simple powerlaw is not 
sufficient at ultra-short periods. This is because for mid- and 
late M-dwarfs (J-K) colour is nearly constant with spectral 



type and as a result the period-colour relation breaks down. 
We used Baraffe (1998) models to convert the maximum 
binary mass for a certain period to the corresponding in- 
frared colours for a single star with the same Tgff as the 
binary. We plot this for equal-mass binaries as a dashed line 
in the right-hand panel of Figure [101 It is consistent with 
the Deb & Singh relation down to ^0.20 days within their 
la uncertainties, but is a better fit to the data at shorter pe- 
riods. Riaz et al. (2008) performed a spectroscopic study of 
1080 nearby active M-dwarfs and obtained 2MASS colours. 
The grey shaded region indicates (J-K) colours for their M- 
dwarfs classified as M4, showing that the infrared colours 
for the ultra-short period binaries are consistent with obser- 
vations of active single mid-type M-dwarfs. 

6.3 Constraints to binary evolution scenarios 

Our results imply that for M-dwarf binaries the orbital peri- 
ods can extend significantly beyond the proposed period cut- 
off of ^^0.22 days. For our (near-) contact M-dwarf system 
19b-3-06008, the fitted Teff -3150 K, which corresponds to 
a mass of O.15-O.25A10 , implies that there is insufficient an- 
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Figure 10. Left panel: Period histogram of the WTS short period binaries (grey and black filled bars), combined with sources from 
the literature (white bars). Right panel: The black stars indicate our WTS eclipsing binary sample with P <0.23 days, whereas open 
boxes are sources from the literature. The period-colour (J-K) power-law relationship for contact binaries from Deb & Singh (2010) is 
indicated by the solid black curve. The dotted curves show the 1 and 3a boundaries on the power-law index of the correlation. The 
dashed curve shows the theoretically predicted colours for periods at Roche lobe overflow using Baraffe IGyr single star models. The 
grey band indicates the range of (J-K) colours for a sample of nearby active M-dwarfs classified as M4 (Riaz et al. 2006). 
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Table 2. Source property table for our red short period sample, P <0.30 days and (r-i)>0.5 and (i-z)>0.25, and our blue sample with 
P <0.23 days. The epoch column shows MJD-245000. The ratio d^ijdx indicates the secondary eclipse depth over the primary eclipse 
depth in the phasefolded J band lightcurve, determined by median averaging the data in a phase range of ±0.01 around phase 1.0 and 
0.5. Quoted magnitudes are SDSS DR7 ugriz and WFCAM JHK, rounded off to two decimals. Our full SEDs have 8 band grizYJHK 
photometry. The state column indicates the estimate from Fourier fitting (Section 4.1) of the morphological state of the binary ('D' for 
detached, 'SD' for semi-detached and 'C for contact.). 



gular momentum loss over the Hubble time within the AML is rather near-equal-mass, unless the secondary is extremely 

framework to explain the observed orbital period with start- bloated in radius with respect to its mass. Also, a mass-ratio 

ing periods of Po=2.0-1.5 days. For the Po=1.0 day model of 0.3 and a primary mass of O.2A^0 would bring the sec- 

and a mass-ratio of ^0.3 the contact timescale is of the or- ondary mass uncomfortably close to the hydrogen burning 

der of the Hubble time. However, the near-equal ratio of sec- limit, whereas an intrinsically dark and low-mass secondary 

ondary to primary eclipse depths suggests that this binary such as a brown dwarf would be unlikely to produce the ob- 
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served secondary eclipse depth and the out-of-echpse varia- 
tions. Again, for the other three ultra-short period binaries 
only the Po = 1.0 day model can reproduce the observed or- 
bits for the estimated primary mass, under the requirement 
that we assume Mi > 0.35Mq and/or an extreme mass ra- 
tio, the latter which may be the case for 19h-3- 14922 (where 
we estimate that q could be as low as 0.5). Given our esti- 
mates of the minimum birth separation for M-dwarfs (Figure 
[2|) , it seems unlikely that these binaries were formed at such 
short periods. 

Another possibility is that the evolution of M-dwarf bi- 
naries is faster than expected. M-dwarfs are known to be 
active and flaring stars. West et al. (2011) find an activity 
fraction of 40-80% for M4-M9 dwarfs, which implies that 
scaling solar type stars to M-dwarfs is not trivial. It is pos- 
sible that because of the (near-)convective nature of mid- 
to-late type M-dwarfs that the ultra-short periods of our 
binaries, when synchronised with the stellar rotation, cause 
the magnetic field lines to be significantly twisted (the 
dynamo; Radler et al. 1990). In other words, the topology 
of the magnetic fields could be significantly different from 
that of solar- type stars, which may directly affect the over- 
all activity and as a consequence the rate at which angular 
momentum is lost. 

Alternatively, the formation mechanism for M-dwarfs 
may be different from that of earlier type main sequence 
stars. It is possible that during the pre-main sequence phase 
an excess amount of angular momentum is removed which 
accelerates the orbital evolution. Possible sources of such en- 
hanced evolution could be dynamical interactions with other 
stars, accretion of the surrounding material and/or interac- 
tion with a circumbinary disk (Pringle 1991; Artymowicz et 
al. 1991; Bate & Bonnell 1997). The hydrodynamical cal- 
culations of binary formation by fragmentation of Bate et 
al. (2002) indicate that proto-binaries potentially form with 
large separations (^10 AU) and go through a phase of accre- 
tion and orbital evolution towards tighter orbits 1 AU). 
Less is known about the later phases of orbital evolution. In 
this paper, we have demonstrated that periods of ^1.0 day 
would be required at early age to explain the observed ultra- 
short periods for our M-dwarf binaries, if magnetic braking 
is the dominant mechanism by which low-mass stars evolve 
over Gyrs timescales. 

Recall that Jiang et al. (2011) concluded that contact 
M-dwarf binaries (primary mass ^O.GSM©) are short-lived 
and thus detections of them during this 'special' phase of 
their evolution would be rare. This is because they predict 
an instability in mass-transfer at the point at which the pri- 
mary fills its Roche lobe which merges the system almost in- 
stantaneously. In combination with the findings by Ge et al. 
(2010), who show that the instability of mass transfer occurs 
promptly upon Roche lobe filling of the primary (a semi- 
detached binary) for stars with a significant convection zone, 
it seems that the existance of our new near-contact M-dwarf 
systems, 19b-3-06008 and 07g-3-05744, among -10,000 M- 
dwarf sources in the WTS Survey, is unusual. However, de- 
tached short-period low-mass binaries are allowed in this 
model, because the angular momentum loss timescale is 
significantly shorter than in the Stepien AML model for 
M-dwarfs. Our detached ultra-short period systems, 17d-3- 
02440 and 19h-3- 14922, thus appear to fit within the Jiang 
et al. model predictions. 



7 CONCLUSION 

In this paper we have presented the results of an exten- 
sive search for short-period low-mass eclipsing binary sys- 
tems in the near-infrared J-band lightcurves of the WF- 
CAM Transit Survey. Probing over 10,000 M-dwarfs down 
to J=18 we report on the discovery of four ultra-short pe- 
riod (P^0.18 days) M-dwarf binaries. Their periods are sig- 
nificantly shorter than of any other known main-sequence 
binary system. All four are below the observed sharp period 
cut-off at P — 0.22 days as seen in binaries of earlier type 
stars. The shortest-period binary consists of two M4 type 
stars in a P = 0.112 day orbit. In total, we find 25 binaries 
with orbital periods shorter than 0.23 days of which 9 are 
identified as likely M-dwarf systems through their broad- 
band colours or spectra. Of these 9 systems, 5 are indicated 
as potential contact systems. These detections pose a direct 
challenge to popular theories that explain the evolution of 
short period binaries by loss of angular momentum through 
magnetised winds, or by unstable mass- transfer, which pre- 
dict timescales that are either too long for their formation, or 
timescales that are too short to observe them in the contact 
phase. Our discovery of binaries with significantly shorter 
orbital periods than the previously observed cut-off implies 
that either these timescales have been overestimated for M- 
dwarfs, e.g. due to a higher effective magnetic activity, or 
that the mechanism for forming these tight M-dwarf bina- 
ries is different from that of earlier type main-sequence stars. 
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Figure 11. The red sample of candidate eclipsing binaries, which satisfies the colour constraints (r-i)^0.5 and (i-z)^0.25, ordered by 
increasing period P ^0.23 days. We plot two full binary orbits. The boxes in the lower right corner of each plot show a three colour 
SDSS composite image of each source. For 17d-3-02440, which does not have SDSS coverage, we show the WFCAM J band image. A 
colour version of the thumbnails is available in the online paper. 
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Figure 12. Red sample with P ^0.23 days continued. 
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Figure 13. Red sample with 0.23< P ^0.3 days. 
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Figure 14. Red sample with 0.23< P ^0.3 days continued. 
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Figure 15. Blue sample with P ^0.23 days. 
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03e_4_2972 period = 0.2 1 931 536 days epoch= 1 1 55.5429 03e_4_3812 period = 0. 22 1 89564 days epoch = 534. 24769 
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19h_1_10160 period = 0. 22338883 days epoch = 374.31 584 1 9c_3_1 2555 period = 0. 2252 1 823 days epoch = 324. 53337 
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Figure 16. Blue sample with P ^0.23 days continued 



Four ultra-short period eclipsing M- dwarf binaries in the WFCAM Transit Survey. 21 
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Figure 17. Blue sample with P ^0.23 days continued 
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Table 3. Table showing the sky coordinates for our WTS short-period echpsing M-dwarf binaries: we hst right ascension (RA) and 
dechnation (DEC) in the format /ims(hours, minutes, seconds) and dms(degrees, minutes, seconds). 



